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Using density-functional theory we compute the energy barriers and attempt frequencies for surface 
diffusion of Ag on Ag (111) with different lattice constants, and on an Ag adsorbate monolayer on 
Pt (111). We find that the attempt frequency is of the order of 1 THz for all the systems studied. 
This is in contrast to the so-called compensation effect, and to recent experimental studies. Our 
analysis suggests that the applicability of simple (commonly used) scaling laws for the determination 
of diffusion and growth parameters is often not valid. 
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The determination of the microscopic parameters that 
govern the evolution and morphology of a growing struc- 
ture during epitaxy is the focus of a large number of 
theoretical and experimental studies. However, exper- 
imentally deduced values for the diffusion barriers and 
attempt frequencies often are not measured directly, but 
result from an analysis of experimental data on the basis 
of some theoretical model. For example, it is qualitatively 
well established that the density of islands N ls in growth 
experiments follows the scaling lawtl 
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Here F is the deposition flux, D is the diffusion constant, 
and Cq is a dimensionless quantity that is determined 
by geometrical factors and parameters that describe the 
capture efficiency of an island. Recently, Eq. (|l|) has 
been applied for a quantitative determination of diffusion 
constants D for different, systems from experimentally 
measured island densitiesoa. 

Relation (Q) is derived under several assumptions, par- 
ticularly that adatoms are the only species that can dif- 
fuse over the surface, i.e., dimers and larger clusters 
are immobile, and that there is only one mechanism 
for atomic diffusion, typically assumed to be the hop 
to a nearest-neighbor lattice site. Then the diffusivity 
D is related to the hopping rate of single adatoms by 
D = Do exp(EdfkT), where is the energy barrier 
for surface diffusion and D = r ^ 2 /(2a) is the pre- 
exponential factor. Here Tq is called the attempt fre- 
quency, I is the jump length of the diffusing adatom, and 
a reflects the dimensionality and symmetry of the motion 
(for a square lattice a = 2). It is commonly assumed that 
To ~ 10 12 — 10 13 s _1 , which is a typical surface phonon 
frequency. We also note that Eq. (|l|) is based on the 
assumption that the so-called critical nucleus is i* = 1, 
which means that dimers and small clusters are stable 
against break-up. This is the case for typical growth ex- 
periments as long as the temperaiure is low enough or 
the deposition rate is high enoughrJ. 

In a recent scanning tunneling microscopy (STM) 
study Brune et ala found that for growth of Ag on a Pt 



(111) substrate the size and density of Ag islands changes 
significantly when one monolayer (ML) of Ag (that grows 
pseudomorphically on top of the Pt substrate) has been 
deposited. The density of islands in both cases is differ- 
ent again from the density of Aa islands on top of clean 
Ag (111). The authors of Ref.Q analyzed their experi- 
mental data with Eq. (Q) with the assumption that the 
mechanism for surface diffusion is hopping of a single 
adatom to a nearest neighbor site. They obtained the 
energy barriers and attempt frequencies quoted in Table 
[jj which shows that the attempt frequency found for Ag 
on 1 ML Ag on Pt(lll) is much smaller than that for 
Ag on Ag (111). A lower attempt frequency for a system 
with a lower diffusion barrier is not an unexpected re- 
sult, and is known in the literature as the compensation 
cffcctQ. It can be understood with simple physical argu- 
ments that are illustrated in Fig. |l|. Consider the solid 
and the dashed line that are the potential energy curves 
for a system with a low and with a high barrier. It is ev- 
ident that the curvature of the potential with the higher 
barrier is larger at the bottom of the potential well. An 
atom in a well with a higher curvature vibrates faster, 
and thus the attempt frequency is Jarger. 

In a molecular dynamics studyQ that employed em- 
pirical potentials and the embedded-atom theory it was 
found that indeed for homoepitaxial diffusion on different 
metal surfaces there is a systematic compensation effect, 
i.e., systems that have a larger energy barrier for surface 

TABLE I. Results for Ed and To . All theoretical values are 
with 3^ = 15, except for Ag on Ag/Pt (111), where 3N = 27. 
We considered an Ag(lll) substrate with tensile (t), compres- 
sive (c), and without strain. The substrate Ag(lll) (c) was 
compressed to the lattice constant of Pt. 



substrate 


experiment^ 


theory (this work) 




E d (meV) 


To (THz) 


E d (meV) 


r (THz) 


Ag (111) (t) 






106 


0.25 


Ag (111) 


100 


0.2 


82 


0.82 


Ag (111) (c) 






60 


1.3 


Ag/Pt (111) 


60 


0.001 


60 


5.9 


Pt (111) 


160 


10 


150 
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FIG. 1. Schematic illustration of the compensation effect 
and the anti-compensation effect. When the barrier increases 
(dashed line compared to solid line), the potential becomes 
stiffer, and the vibrational frequencies increase. However, 
when an increase of the barrier is due to an increase in the 
lattice constant, the curvature of the potential might very well 
decrease (dashed-dotted line), and as a result the vibrational 
frequencies decrease. In this simplified picture it is assumed 
that all the vibrational frequencies at the saddle point increase 
(decrease) similarly to the ones at the lattice site, and thus all 
the increases (decreases) of the frequencies in Eq. (1) cancel, 
except one that corresponds to a vibration in the lattice site. 

diffusion have a larger attempt frequency. Values for Tq 
are obtained that vary from 3 x 10 12 s _1 to 5 x 10 13 
s" 1 , while for the same systems the energy barriers for 
diffusion increase from 0.05 eV to 0.7 eV. The magni- 
tude of the effect for To is in the range of typical surface 
phonon frequencies, and it does not compensate com- 
pletely (for typical growth temperatures) the increased 
diffusion barrier. In particular, we note that the attempt 
frequency is always much larger than the reported value 
of T = 10 9 s" 1 for Ag on 1ML Ag on Pt(lll), which is 
about three orders of magnitude smaller than what one 
would expect. Despite this puzzle, we note that previous 
theoretical studies seem to suggest that the system 1 ML 
Ag on Pt (111) behaves rather normal and that, with 
respect to surface diffusion, it behaves essentially like a 
laterally compressed Ag (111) substrateErtl 

The diffusion constant is determined by the diffusion 
barrier and pre-exponential factor, but for the latter no 
ab initio calculations have been reported so far. This 
applies in general, not just to the system of concern in 
this Letter. In view of the above mentioned puzzle, and 
similar enigma that seem to exist for other systems, we 
performed density functional theory calculations of the 
attempt frequencies for Ag diffusion on Ag (111) with 
different lattice constants, and on 1 ML adsorbed on Pt 
(111). We find that in all cases the attempt frequency 
is of the order of 1 THz, so that claims of a large vari- 
ation of Tq are not supported. We then discuss various 
mechanisms which become operative in systems with low 
energy barriers. 



TABLE II. The calculated attempt frequency To for Ag on 
Ag(lll) for different numbers of degree of freedom considered. 



3iV 


3 


15 


99 


To (THz) 


1.55 


0.82 


0.71 



We employ density-functional theory with the local- 
density approximation!!. More details are given in Ref.tl. 
Unless noted otherwise, we use a supercell with 4 layers, 
a (2 x 2) surface cell, and the adatom placed on only one 
side of the slab. Tests have shown that such a cell is large 
enoughS. 

The energetically most favored adsorption sites are the 
three-fold coordinated fee and hep sites, and an adatom 
that diffuses from one site to another will pass over the 
bridge site as the transition state£L3. In principle, the 
effective diffusion constant on a surface with several ad- 
sorption sites within one unit cell needs to be calculated 
as a combination of the different rates for the motion be- 
tween the different sites. However, for all the systems 
studied here the fee site is always slightly favored over 
the hep site (by 5-10 meV), but the barrier to hop from 
an hep site to an fee site is almost identical to the one 
to hop from an fee site to an hep site. This provides 
justification that only the fee site is considered as the 
(stable) adsorption site. Transition-state theory (TST) 
within the harmonic approximation-is used, and the at- 
tempt frequency is calculated fromliil 



The Vj (v*) are the normal-mode frequencies of the sys- 
tem with the adatom at the fee site (bridge site), and 3N 
is the number of degrees of freedom of the system. 

We have calculated these normal mode frequencies by 
first evaluating and then diagonalizing the force-constant 
matrix. The elements of this matrix are obtained by 
displacing each individual atom in the cell and calcu- 
lating the forces that act on all the atoms. An impor- 
tant issue is how large 3N should be for constructing the 
force-constant matrix. The results obtained for differ- 
ent 3N are summarized in Table || for Ag on Ag (111). 
If we change 3N from 3 to 15, i.e., take into account 
the adatom and all atoms in the first layer, the attempt 
frequency changes from 1.55 THz to 0.82 THz. We also 
developed a scheme to include a larger number of degrees 
of freedom, and given in Table || is a value for 3N = 99 
which corresponds to a (4 x 4) surface cell with vibra- 
tions of the adatoms and two layers of the substrate. The 
(99 x 99) dynamical matrix was constructed as follows: 
We use the same force constants as in the 37V = 15 situ- 
ation in the corresponding subspace. For the other force 
constants we use bulk values (for atoms in the second 
and deeper layer), and clean-surface values (for farther 
surface atoms). Tests with a (3 x 3) cell revealed that 
force constants of surface atoms due to the displacement 
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of atoms that are further away than nearest neighbors are 
negligible. Similarly, force constants of atoms in the sec- 
ond layer (and below) that result from the displacement 
of atoms that are further away than its nearest neigh- 
bors can be neglected. We estimate that the influence of 
these negligible force constants on the attempt frequency 
is less than a factor two, and therefore set them to zero. 
The attempt frequency for 3N = 99 decreases very lit- 
tle compared to the 3N = 15 treatment, namely to 0.71 
THz, and thus our result obtained when only the adatom 
is allowed to vibrate-,is correct within a factor of two. 

As argued earlieroB, the main difference between an 
Ag atom on a clean Ag (111) surface and on 1 ML Ag on 
Pt (111) is due to strain in the Ag adlayer. We therefore 
studied the attempt frequencies for Ag on Ag (111) as a 
function of strain. The results are given in Table H. An 
increasing lattice constant implies that the diffusion bar- 
rier increasesEl, so, according to the compensation effect, 
the attempt frequency should increase. Surprisingly, the 
opposite is the case: When the barrier increases from 60 
meV to 106 meV, the attempt frequency decreases from 
1.3 THz to 0.25 THz. We understand this "wrong trend" 
as follows: When the diffusion barrier increases because 
the lattice constant increases, the curvature of the poten- 
tial decreases, so that the attempt frequency decreases. 
This anti-compensating effect is sketched schematically 
in Fig. |l| by the dashed-dotted line. 

Our result for Ag on unstrained Ag (111) differs by a 
factor of ~ 4 from that obtained in the analysis of the 
experimental data. This difference is partially due to the 
fact that the analysis of the experimental data ignored 
that an atom that hops from an fee site to a neighbor fee 
site first visits an hep site. When the difference in adsorp- 
tion energies between the fee and hep sites is very small, 
the adatom needs to make (on average) three times as 
many hopping attempts to get to one of the six adjacent 
fee sites than it needs to hop to one of the three adjacent 
hep sites. Thus, the attempt frequencies given in Ref. ^, 
which are also noted in Table |[ should be multiplied by 
a factor of three. 

For the diffusion of an Ag adatom on 1 ML Ag ad- 
layer on Pt (111) the experimentally deduced attempt 
frequency is more than two orders of magnitude smaller 
than that for the Ag (111) surface, but such a reduction 
is not found in our calculations (compare Table ||) . When 
only the adatom and the Ag layer are included we obtain 
To ~ 3.8 THz (3iV = 15), and upon inclusion of the top 
Pt layer we obtain T ~ 5.9 THz (3N = 27). As men- 
tioned earlier, with respect to diffusion energy barriers, 
the Ag monolayer on Pt (111) behaves very much like a 
compressed Ag (111) surface. We now find that this also 
holds for the attempt frequency. The disagreement with 
the experimentally deduced value, that is three order of 
magnitude smaller, is significant, and we consider this 
disagreement a most important result. It demonstrates 
that the understanding of surface diffusion and growth 
is far from being complete, although often this is not 
appreciated. In this regard we note that Eq. (|l|) is de- 



rived under several assumptions, and that often some of 
them are not met. We will now discuss six different prob- 
lems which may have troubled the experimental analysis 
and/or our theoretical results. 

i) : It is a concern whether the harmonic approxima- 
tion we used is valid for all the systems under considera- 
tion. From the calculation of the forces we find that for 
the systems studied here the harmonic approximation is 
indeed a rather good approximation for displacements up 
to 0.2 A. 

ii) : It is known that TST is only working when the 
following conditions are met: The adatom spends enough 
time at the adsorption site so that it can equilibrate 
there. This aspect is connected to the phononic friction, 
i.e., the strength of the coupling of the adatom vibration 
to the excitations of the substrate. Additionally, once an 
atom has crossed the transition site there should be no 
recrossing before the atom has been equilibrated again. 
While it is clear that the friction should not be too small 
or too large, .a. quantitative criterion when TST is valid 
is not knowno. We expect that equilibrium might not 
be reached when the temperature is too high, but the 
experiments of RefBI where carried out at temperatures 
that are smaller than half of the Debye temperature. 

in) : Diffusion of small clusters on a metal surface can 
be significanto, but Eq. (0) is only valid when the oniy. 
mobile species are adatoms. In fact, it has been shown t3 
that the scaling exponent in Eq. (gj) changes when small 
clusters are allowed to_move and that the island density 
decreases dramatically^. This effect cannot explain the 
discrepancy we find for the diffusion on Ag/Pt (111): If 
cluster diffusion were important on Ag/Pt(lll) and not 
(or less) on clean Ag (111), the pre-exponential factor for 
diffusion on Ag/Pt(lll) should be bigger, not smaller, 
as observed in the experiments. Still, we do not rule out 
that cluster diffusion plays a role, because it is possible 
and likely that the diffusivity of cluster diffusion is (very) 
different to that of single-atom diffusion, and we also note 
that the exponent in Eq. ([l]) may change. We are not 
aware of the direction and size of this effect, but note 
that cluster diffusion may be particularly important for 
strained systems and thus Ag diffusion on Ag/Pt(lll) is 
indeed a good candidate. 

iv) : It has been shown that on metal surfaces the dif- 
fusion-length may be a multiple of the lateral lattice con- 
stantOlij. This becomes more relevant when the corru- 
gation of the potential energy surface jsdess pronounced 
and the diffusion barrier becomes smallEj, as it is the case 
for Ag on Ag/Pt(lll) and Ag on a compressed Ag (111) 
substrate. The adatom becomes more delocalized and we 
speculate that rather than hopping over the surface from 
one lattice site to an adjacent one, the adatom is "flying" 
over the surface over several lattice sites. For this process 
TST is strictly no longer applicable. Moreover, it is con- 
ceivable that Eq. ([!]), particularly its exponent, changes 
when several diffusion mechanism are operative. We are 
not aware of the size or direction of such an effect, but 
under the speculative assumption that the exponent is 
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reduced by 20%, the discrepancy between the results of 
Refil and our calculations could be explained. The anal- 
ysis of the experimentally measured island density would 
then yield an activation barrier, that it 20% larger, which 
is supported by experimentall3 and theoreticaO results 
for the double-jump of a Pt adatom on Pt(110). This 
would also explain the good agreement between our cal- 
culated and the experimentally deduced diffusion barri- 
ers. Additionally, an exponent that it 20% smaller yields 
an attempt frequency of 1 THz, in agreement with our 
calculated value. 

v) : We also note the possibility of diffusion by atomic 
exchange over long distances. The Ag adlayer is com- 
pressed and we speculate that atomic exchange could 
work such that an atom of the Ag adlayer pops up, be- 
coming an adsorbate atom. Then the vacancy moves to- 
wards another adatom (e.g. a long row of Ag adatoms 
in the Ag adlayer may move). Similar effects have been 
observed in embedded atom simulationsEj. 

vi) : Another problem is knowledge of the dimension- 
less quantity Co in eq. (]!]). An uncertainty in Co of one 
order of magnitude translates into an uncertainty for T 
of three orders of magnitude when i* = 1. The value of 
Co is essentially determined by the surface geometry and 
capture efficiency of islands, and with the self consistent 
approach of Bales and Chrzannil one obtains Co — 0.25, 
which has been used in RefJ-3. A recent STM study for 
Pt on Pt(UJ) used an alternative method to determine 
i?d and r E3. It is found that Jlp = 5 THz (an expected 
result), and iom the data inE3 we deduce Co — 0.28. 
Earlier workEj suggested that Co should be of the or- 
der of unity. It is known that often there is funneling 
of adatoms toward step edges. This attraction will be- 
come more significant for systems with low diffusion bar- 
riers [as it is the case for Ag on 1 ML Ag on Pt (111)], 
and when surface states are important and affected by 
adatoms and/or steps. 

We conclude that attempt frequencies of single-atom 
hopping diffusion are indeed of the order of 1 THz and 
reported deviations from this value, as. for example the 
reduction by three orders of magnitudeo, are only appar- 
ent and actuated by the low energy barrier of diffusion. 
This implies that adatoms should be described by a more 
delocalized nuclear wavefunction. Long-range adatom- 
adatom and atom-step interactions become noticeable, 
and atomic diffusion proceeds more by "flying" over sev- 
eral lattice constants rather than by hopping to nearest 
neighbor sites. For strained systems also cluster diffu- 
sion may be enhanced. The importance of these effects 
is difficult to assess by theory, but we trust that more 
experiments will help to clarify the situation. 
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